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We present preliminary results of a new approach to the study of the pion-pion system in the 1=1 
channel. The Bethe-Salpeter wave function of the two-pion system is computed on the ground 
state and the first excited state. From these, we attempt to extract an interaction kernel (potential) 
which can then be used to extract observables such as the phase shifts. In a first trial, we use 
rather large pion masses m K ~ 1 .05 GeV and m n ~ 0.68 GeV which do not allow rho decay. 
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1. Introduction 

In recent years, the 1=1 two-pion system has attracted a lot of attention in lattice QCD. The 
increase in computational power has finally allowed to generate fully dynamical QCD gauge con- 
figurations at quark masses low enough for the rho resonance to be observed, promising a better 
understanding of complex hadron processes from first principles. Recent studies [[]], |2[ B[ all use 
the de facto standard method, which is to apply Liischer's formula [Q] or its extension to moving 
frames [H], in order to relate the finite-size energy spectrum to the infinite volume phase shifts. The 
main difficulty faced by this method is that one can only extract the phase shifts at a few energies 
on the lattice, making it difficult to reconstruct the continuous energy range and thus the physical 
parameters of the system, especially around the resonance. 

An alternative method has been recently introduced for the study of the nucleon-nucleon sys- 
tem |7|] then successfully extended, in particular, to various baryon-baryon systems []8|]. The method 
relies on the fact that the phase shifts can be obtained from the asymptotic behaviour of the Bethe- 
Salpeter (BS) wave functions. An effective energy-independent, non-local potential can then be 
introduced to account for the energy dependence of the BS wave functions, and as a result, the 
energy dependence of the phase shifts. 

This paper reports on our first attempt to apply this potential method to the 1=1 two-pion 
system. The meson masses considered here do not allow for the rho meson to decay, the goal being 
to test the viability of the method in this channel before applying it to the study of the resonance. 

2. BS wave function and potential 

To describe two pions in the isospin 1 = 1 channel, we use the following operator 

7T7r( P ) = -= [7r(p)7T + (-p) - 7r+(p)7r(- P )] , (2.1) 

where ^ ± are local interpolating operators for the pions. The Bethe-Salpeter (BS) wave function 
in the center of mass frame is then defined in this channel as 

%(r) = /(03^ rp (O|^(p)|n) (2.2) 

with \n) an eigenstate of QCD with the required quantum numbers. The BS wave function presents 
an asymptotic behaviour [^] which allows to extract the scattering phase shift at the energy of the 
eigenstate, in exactly the same way as the wave function in quantum mechanics. 

By inversion of the energy-dependence of the wave functions [gj, one can define a non-local 
energy-independent potential U such that the wave functions satisfy, for all eigenstates \n) with 
energies below the inelastic threshold, the Schroedinger-like equation 

(V 2 + kl)y n (r)= mjl j d 3 r' U(r,r')y n (r'). (2.3) 
where the energy of the eigenstate \n) is E n = 2 J + m\. 
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Since we only have access to a limited number of such BS wave functions on the lattice, we 
can only obtain an approximation of the potential, such as a truncation of its derivative expansion. 
The velocity expansion of the non-local potential in this channel is 

U(r,r') = [Vo(r)+V 2 (r)V 2 rl + . . .]8(r' - r). (2.4) 

From the approximate potential, we can obtain approximate BS wave functions at all energies 
below the inelastic threshold and extract phase shifts from them. The accuracy of the results thus 
depend on the convergence of the derivative expansion. 



3. Wave functions on the lattice 

3.1 Eigenstates and correlation matrix 

We note G(t) a functional of the quark fields at time t, and 6 its corresponding operator in the 
Heisenberg picture. In the limit of an infinite lattice in the time direction, we have 

(<7i(*Wb)> =£(0|^i|«)(n|^ 2 |0)^ £ "^. (3.1) 

n 

where the brackets in the left-hand side denote the expectation value in lattice QCD. 

Let be a set of N linearly-independent operators and M n j = (n\ ffj\0) the N x N matrix of 
their mixing with the N lowest eigenstates. Their correlation matrix G is defined as 

Gij{t,t ) = (^(tWj(t )) =M^D(t-t )M + O(e- E ^-^) (3.2) 

where D is diagonal with components D nn (t — ?o) = e~ E "('~' \ 

Diagonalizing G~ l (t' ,to)G(t,to) for several (?,?') pairs yields the energies, by fitting of the 
eigenvalues with the form e E ' , ( t '~'\ and the inverse mixing as the eigenvector matrix (up to 
a normalization of the columns). This requires t and t' to be sufficiently separated from to so that 



contributions from eigenstates higher than \N) vanish in (3.2). 

The determination of M~ l and the energies allows to compute, for any operator & and ei gen- 
state n <N, the matrix element 

(0\^\n}=e E "^^(M-% 1 (ff(t)W i (t )) (3.3) 

i<N 



and thus the BS wave functions as seen in the definition (|2.2|). 
3.2 Source operators 

The pion-pion 1=1 channel contains both pion-pion scattering states and the rho meson. In 
the pion mass region we investigate {m % = 1.05 and 0.68 GeV), we expect the ground state to be 
the rho meson, the first excited state to be the pion-pion scattering state with the lowest non-zero 
momentum allowed on the lattice, and other eigenstates to have energies large enough to only 
consider the first two (N = 2). 

To approximate the pion-pion state we use the operator G\ = nn(p) with momentum p = 
(L being the spatial extent of the lattice) and for the rho meson the operator 

02 = p = -^Y [u(x)a-yu(x) - d(x)a-yd(x)] (3.4) 
v2 x 

with a polarization taken parallel to that of the relative momentum of the pions, a = e3. 
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(jwr(q;f)p(*b)> = 




- (q <-> -q) 



Figure 1 : Decomposition in Wick contractions of the correlation functions corresponding to kk — > 
Tin and p — > Tin, appearing both in the correlation matrix (with q = p) and the wave function. Time 
goes upward. 




(a) (b) (c) 



Figure 2: Computation method for some contractions. Springs link pairs of points which are pro- 
jected one on the other by summing over stochastic noises. Open circles are explicit summations. 
Straight (resp. broken) arrows are direct (sequential) propagators. 



3.3 Correlation functions 

The correlation matrix and the wave functions require the evaluation of correlation functions 
of the form of (3T ), some of which are illustrated in Fig. [T| For the wave functions, the substitution 



q f-> — q translates to parity conjugation, so we only need to compute the first parts. The diagrams 
are computed, following [jl|], using adequate contractions of direct and sequential propagators using 
stochastic noises £j as sources 

e(x, ? |q, ?5 ,^) = I^ 1 (x,?;y,^)[^ q - y ^(y)] (3.5) 

y 

W(x,f|k,f 1 |q,f 5 ,^)=££>- 1 (x,?;z,f 1 )[e' k - z e(z,fi|q,f5,^)]. (3.6) 

z 

Fig. ^| shows an example. The propagators are contracted explicitely at the sink (upper part 
in the diagram) and implicitely at the source (lower part) by noise projection. In Fig. [I] appear 
rectangle- and triangle-like diagrams. They can be computed in the same way as the previous ex- 
ample using sequential propagators, cf. Fig. |b|. However, while the momentum can be introduced 
freely at the explicit summation (empty circle), sequential propagators have definite intermediate 
momenta. This means that for the wave functions, which requires all possible sink momenta q, 
we need to compute as many sequential propagators. To remedy this, we introduce for the wave 
functions another stochastic noise at the sink, cf. Fig. which allows to choose the two momenta 
at the sink independently of the computation of the propagators. 
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Figure 3: Wave functions for the first (left) and second (right) diagrams of the nn — > Tin correlation 
function in Fig. [j] (upper). Normalized such that the total wave function has a norm 1. Computed 
at t-t = 12 on the set K ud = 0.1347. 



4. Numerical setup 

The preliminary results presented here were computed using the Nj ■ = 2 + 1 full QCD gauge 
configurations of ILDG/JLDG generated by the CP-PACS and JLQCD Collaborations on a 
28 3 x 56 lattice with a RG improved gauge action at j3 = 2.05 and a 0(a) improved Wilson quark 
action with csw = 1.628. The lattice spacing is a = 0.0685 fm which makes for a lowest non-zero 
momentum of p = 2%jL = 0.65 GeV. We compare our results on two sets of configurations with 
light quark hopping parameters K ut i = 0.1347 and K u d = 0.1356, keeping K s = 0.1351 fixed. 

The configurations with K u d = 0.1347 have meson masses m n = 1.05 GeV and m p = 1.37 
GeV. Those with K u d = 0.1356 have masses m % = 0.68 GeV and m p = 1.10 GeV. In both cases, 
the lowest energy of two free pions in the center of mass frame, 2\Jm % + (2n/L) 2 , is significantly 
larger than that of the rho meson at rest. 

The quark propagators are computed with temporal Dirichlet boundary condition. We use 
U(l) stochastic noises, 6 at the source and 20 at the sink. Wave functions are projected in the 
representation of the cubic group. Statistical errors are computed using the jackknife technique 
although 2-dimensional plots are shown without error bars for clarity. 



5. Preliminary results 

We have seen that the wave functions ( j2~2| ) on the eigenstates are obtained as combinations 
(33) of the wave functions computed with the interpolating source operators 7in(p) and p, them- 
selves computed as sum of Wick contractions ("diagrams"). The combinations are obtained by 
diagonalization of the correlation matrix. 

Figure || shows the contribution of the two kind of diagrams appearing in the %% — > %% wave 
function. The left one, corresponding to the "parallel" diagram, is close to the free wave function. 
The right one, corresponding to the rectangle diagram, exhibits a very peaked and short-ranged 
behaviour. The Uiangle diagram, Fig. [T] (lower), has a wave function very similar (up to a normal- 
ization) to that of the rectangle one. The rho meson being the ground state, the quark-antiquark 
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pair propagating in the rectangle and triangle diagrams from t to to can be thought as forming a rho 
meson, which could explain why the two diagrams' wave functions are similar and short-ranged. 

The ground state's wave function, Fig. |] (left), can be obtained using either source operator 
by saturation at large enough time separation. Taking 7T7r(p) as source operator, we see that the 
dominant contribution as time separation increases is from the rectangle diagram. 

The first excited state's wave function is shown Fig. |] (right). We see that the dominant 
contribution is this time coming from the parallel diagram. The wave functions is obtained with a 
linear combination of the two source operators, which has for effect the cancellation of the peaked 
short-range contribution between the rectangle and triangle diagrams. However, while the signal 
from the ground state wave function is cancelled, the statistical noise remains and grows as e^i'-'o) 
with AE the energy difference between the two lowest eigenstates. 

An approximate potential is obtained by inverting the Schrodinger equation (2.3) with the BS 
wave functions computed on the lattice as input. The wave functions in Fig. f| unfortunately do 
not allow such a computation. The ground state wave function (left) is sharply peaked around the 
origin, leading to huge discretization errors when taking finite-difference Laplacian operator. The 
first excited state wave function (right) is extremely noisy due to large energy separation between 
the two lowest eigenstates and the noise is further enhanced by taking the Laplacian. 

Using the fact that the main contribution to the first excited state wave function is from the 
parallel diagram and that the other diagrams should only contribute to the short-range part of the 
potential, we show in Fig. |5| the effective central potential computed using only the parallel dia- 
grams, on both sets of hopping parameters. We see that a simple Yukawa fit is in good agreement 
to the data even at surprisingly short range. The masses in the Yukawa fit are 1.53(9) GeV and 
0.94(17) GeV, with corresponding rho masses of 1.37 GeV and 1.10 GeV, respectively. 



6. Summary and outlook 

We have shown preliminary results of the application of the potential method to the 1=1 pion- 
pion system. The method, which has been successful in the study of baryon-baryon systems, 
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Figure 5: Central potential using only parallel diagrams, t — to = 12. Fit by a Yukawa potential. 



encountered difficulties in this particular setup. First, the ground state being the rho meson, the 
wave function is very short-ranged and the computation of the potential leads to large discretization 
errors. Then, while the first excited state is a scattering state and likely to be well described by a 
potential, it is difficult to extract due to the large energy difference. 

However, approaching the problem from a different perspective, the present results shed a new 
light on the qualitative understanding of the system. Furthermore, the above problems may be 
solved in the region where the rho meson is a resonance and not the ground state, since the scat- 
tering state will be simply extracted by saturation and the short-ranged component should become 
less important. In this case, the potential method could lead to competitive quantitative results. 
Further study at smaller pion masses will confirm or invalidate this expectation. 

Numerical computations in this work were carried out on SRI 6000 at YITP in Kyoto Uni- 



versity. We are also grateful for the authors and maintainers of CPS++ pi]], of which a modified 
version is used for measurement done in this work. 
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